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This talk presents a brief review of some theoretical aspects related to the 
physics of color glass condensate (for detailed reviews see, e.g., [HI2] with the 
main focus on the results obtained in EJ Ej ■ In its essence the color 
glass condensate physics is based on the quasiclassical perturbative approach 
to high energy QCD, in which contributions of all orders of a s ln 1/x, where 
x is a fraction of the hadron (nuclear) longitudinal momentum carried by the 
considered gluon mode, and of all orders in gluon density are resummed. 

The main issue of the perturbative approach to high energy QCD is that 
of uitarity. Resummation of leading energy logarithms 0(a s \nl/x) in the 
linear approximation in gluon density through the BFKL equation for the 
(unintegrated) gluon density 4>(Q]_,x) jH]: 



The growth of the gluon density and the resulting growth of physical cross- 
sections with energy (1/x ~ y^) is powerlike - in clear violation of the 
Froissart bound a ~ 1/ m 2 In 2 1/x following from the requirement of unitarity 
of the theory. 

Let us note, that besides color glass condensate approach there exist two 
other major lines of research addressing the question of unitarization of high 
energy perturbative asymptotics of QCD, in particular: 

• An effective theory including, in addition to quarks and gluons, new 
degrees of freedom - reggeons j2l CD| • 

• Generalization of operator product expansion formalism 

The notion of color glass condensate has been suggested with the following 
two analogies in mind: 

• The small -x gluon modes are characterized by large occupation num- 
bers. Therefore one is tempted to describe physics of gluons in terms 
of classical gluon fields, corresponding to the condensation of gluonic 
modes. 

• Averaging over color charge density is similar to averaging over disorder 
in random magnetics - spin glasses. 




(1) 



Tree level: McLerran-Venugopalan model 
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Foundations of color glass condensate physics lie in the quasiclassical tree- 
level description of the wave function of large nuclei developed by L. Mcler- 
ran and R. Venugopalan The gluonic content of this wavefunction is 

described in terms of fast (having large k + ~ P + , where P + is a longitudinal 
momentum of the nucleus) static sources characterized by the color charge 
density p a (x~,x±) and slow (small k + <C P + ) 

[D U ,F^] = ^ + Pa (x-,x ± ). (2) 

The solution of (J2J) reads: A + = A~ = and 

AHx-,x±) = -U(x- J x ± )d i UUx-,x ± ) , (3) 
9 

f W(x~,x±) = Pexp jig f£- dz~ a(z~,x±)\ 
\ -V^ = WpU 

In what follows we shall also need the following approximate expression for 
W that often enters the expressions arising in calculations of quantum cor- 
rections: 

W(x~, xj_) = Pexp jig J dz~ a(z~ , ~ 6(x~) £tf(x±) + 9(— x~), 

(4) 

Calculation of physical observables is then done through computing the tree 
level correlators on the classical solution A % {x~ ,x±): 

(Ai(x\x)Ai(x\y)...) A = fvp W A [p]Ai(x)Ai(y)--- (5) 



where W\ [p] is a weight functional for the ensemble of sources at the (longi- 
tudinal) scale A. 

The main result of the tree-level approach is that already in this ap- 
proximation the infrared growth of gluon distributions is tamed by non- 
linear corrections |3J EJ. This happens at the new scale of the theory - 
the saturation scale Q s . Let us illustrate the saturation phenomenon by 
writing down the asymptotical expressions for the gluon structure function 
G(x,Q 2 ) = f Q2 dki<j)(x,ki)/kl: 



xG(x,Q 2 )(k 2 ± ^>Q 2 s ) ~ N c ^\n- 
xG(x,Q 2 )(k 2 ± «Q 2 ) ~ ^-^vr^lngl 



A 2 

ly QCD 
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Quantum corrections: renormalization group in r = lnl/x 

The starting point for constructing a theory allowing to compute quan- 
tum corrections to McLerran-Venugopalan model is writing down an effective 
action for high energy QCD in terms of classical sources p and gluon fields 
A 1 . Such an action was first written in [3|: 

S[A,p] = - J d i x~F^Fr + ^ J d 3 £Tr{p(f)VPV-oc(z)} = S YM + S W . 

(6) 

where 



^00,-00 [A~] (x) = T exp 



ig J dx + A (x) 



(7) 



Note that the effective action (jHJ) includes the standard Yang-Mills piece 
Symi as well as a gauge invariant generalization of the abelian eikonal vertex 
/ d A xp a A~. As shown in this action allows to reproduce the BFKL 
equation .8] in the small-density limit. 

The considered effective theory contains two averagins: over the config- 
urations of the classical source p at the scale A with the weight functional 
Wa[p] and usual quantum averaging over the quantum fluctuations over the 
gluon fields. For example, the two-point correlator reads 

« t^o »-Jv, w M { ^ffiX™ } ■ m 

As shown in jUEllZj, the quantum evolution of gluon correlators is com- 
pactly parametrized by the evolution of the weight functional Wa[q] (it turns 
out convenient to make a change of variables from p to a related through the 
solution of the above- described classical equation of motion): 

8W T \a] fl P s .1 

The evolution © includes new virtual and real kernels v and v. Explicit 
computations [HJ Ej lead to the following simple expressions for them: 



v 



(2ir)*(x ± -z ± : 
1 r d 2 z ± (x 1 — z l )(y l — z l ) 

x {l + Q\x ± )Q(y L ) -Q\x ± )Q(z ± ) -Q\z ± )Q(y L )} ab (10) 
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As shown in the coefficients v ah and v a are related to each other by the 
following functional identity 

1 r 2 Su ab {x ± ,y ± ) 

2J d ^ Sa»(y ± ) =vM (U) 

The functional identity (jTTj) allows to rewrite the evolution equation Q in 
the elegant Hamiltonian form [TK] 

0W r [a] f /"d 2 ^ 



d 2 x± z % — x l r +/ . . ,\ 5 



In the weak field limit (weak a) one recovers the BFKL equation 8j. 

Evolving the n - point correlator with the BFKL kernel JCbfkl leads 
to the BKP equations. 



Of special interest is the evolution of V(x±,y±) = Tr \pJ(x±)£l(yj 
It is described by the Balitsky-Kovchegov equation [TTl ITBj : 

dV(x ± ,y ± ) r 2 (x ± - y ± ) 2 

dr 2tt 2 7 Z± (x ± - z ± ) 2 (y ± - z ± ) 2 [ cVxy VxzVzy) 

(13) 

Using Eq. IjlHj) one can write an evolution equation for the scattering 
amplitude of color dipole on the (gluonic field of the) target J\f(x±, y±): 

Af{x ± ,y ± ) = ±-Tr(l-tf(x_ L )n(y ± j) (14) 

m{x ^ y±) = fd 2 z ± - ( {x± " yxf [N xz + N zy - N xy - N xz M zy ] 

Or J (x± - z ± ) 2 (y ± - z ± ) 2 

Note that in the Balitsky-Kovchegov equation (|T3|) the contribution 
quartic in Q is factorized. 

At present the Balitsky - Kovchegov equation (JT3j) has also been derived 



in 



Regge formalism, ^7j 
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— QCD-Regge formalism jTH] 



Some important results 

Let us now mention some important results obtained in the course of 
developing the color glass condensate physics. 

• Nonlinear effects taken into account in (j!3|) ensure unitarity restoration 
for the scattering at fixed impact parameter only [T9| l20| l2Tj . The 
inelastic cross-section is still characterized by the powerlike growth 

&inel = ^Rtarget + 2n R tar getXoe UT (15) 

The origin of the unitarity violation is the powerlike decay of the gluon 
fields generated by the multipole configurations under consideration. 
To enforce the Froissart bound one has to ensure the exponential decay 
of the gluon fields in the impact parameter plane, which is equivalent to 
the existence of the mass gap of the theory. Obviously, no perturbative 
framework can describe such an effect. 

• The functional form of the above- described tree-level solution remains 
unchanged in the course of RG evolution. All changes can conveniently 
be parametrized by evolving the saturation scale Q 2 S [221 123 121] : 

Q 2 s (b±,r) = Q 2 s (b ± ,T )e Xa ^-ro) 
N(x±, y±) = l-exp[-r 2 ± Q 2 s (b ± ,r)] (16) 

Physically this situation can be described as a geometrical scaling of 
observable cross-sections, when all physical quantities depend on energy 
through Q 2 (- ■ ■ , r) in some range of transverse momenta: 

o = o (Q 2 /Q 2 S ) , Q 2 < Qt/A 2 Q cD (17) 

• Phenomenological implications of the above - described "Q 2 " - physics 
are being actively studied for ultrarelativistic heavy ion collisions, deep 
inelastic scattering on nuclei, etc. 

Conclusions 
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Let us conclude with summarizing once again the main points touched in 
the talk and formulate some major unresolved issues: 

• Quasiclassical RG formalism provides a compact and transparent view 
on nonlinear effects in QCD at high energies 

• Nonlinear effects in leading logarithmic approximation provide unitar- 
ity restoration at fixed impact parameter only. This leaves us with a 
number of important unsolved questions. Two of them are: 

— How many pomerons does one need to describe high energy QCD? 

— Does perturbative unitarity exist? 

The work was supported by RFBR grant 04-02-16880 and the Scientific 
Schools Support Grant 1976.2003.02. 
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